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The interlayer magnetoresistance of a quasi-one-dimensional Fermi liquid is considered for the 
■ case of a magnetic field that is rotated within the plane perpendicular to the most- conducting 

f***) ' direction. Within semi- classical transport theory dips in the magnetoresistance occur at integer 

, "magic angles" only when the electronic dispersion parallel to chains is nonlinear. If the field 

direction is fixed at one of the magic angles and the temperature is varied then the resulting variation 
of the scattering rate can lead to a non-monotonic variation of the interlayer magnetoresistance with 
temperature. Although the model considered here gives a good description of some of the properties 
of the Bechgaard salts, (TMTSF) 2 PF 6 for pressures less than 8 kbar and (TMTSF) 2 C10 4 it gives 
I , a poor description of their properties when the field is parallel to the layers and of the intralayer 

transport. 

^ \ PACS numbers: 74.70.Kn, 72.15.Gd, 71.10.Hf, 74.20.Mn 

i: 

In spite of intensive research over the past decade the nature of the metallic state in low-dimensional strongly 
correlated materials is still poorly understood. Widely studied materials include cuprate and organic superconductors 
[0,1). Many of the properties of the cuprates cannot be understood within the Fermi liquid picture that has so 
successfully described conventional metals ||. Although some properties of the quasi-two-dimensional molecular 
crystals, k-(BEDT-TTF) 2 X (Reference (§) and the quasi-one-dimensional Bechgaard salts [Q (TMTSF) 2 X can be 
explained within a Fermi liquid framework others cannot. A particular challenge is understanding the dependence of 
the magnetoresistance of the Bechgaard salts on the direction of the magnetic field, especially (TMTSF) 2 PF6 under 
pressures of about 10 kbar ^j?]]. The different angular-dependent magnetoresistance effects in quasi-one-dimensional 
' metals are known as the Danner magic angle (or Lebed) [p|-^4|, and third angular effects jlfj, depending on 
whether the magnetic field is rotated in the a c, b c, or a b plane, respectively. (The most- and least-conducting 
• directions are the a and c axes, respectively). The magic angle effect is the most poorly understood of these effects 
and is the focus of this paper. If 9 is the angle between the magnetic field and the c axis then at the "magic angles" 
^ • given by 

tan0 = -- ±P,q= 1,2,3,.... (1) 
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I. INTRODUCTION 
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where b and c are the lattice constants in the b and c directions, Lebed predicted dips in the threshold field for 
formation of a field- induced spin-density- wave (9|. Although these dips are not observed jio), features are seen, 
mostly at p/q — 1,2 in the torque JlJ] and in all components of the resistance [(t|[12] pd|Jlrj -[l8|| ■ 

A wide range of physical mechanisms have been proposed to explain these effects including commensurability effects 
, changing the electron-electron scattering rate pl| , semi-classical transport fll9|| , complicated band structures |po| , pl| , 
hot spots on the Fermi surface , cold spots on the Fermi surface |2^j , electron-electron interactions |24|] , non- Fermi 
liquid effects and magnetic field induced changes in effective dimensionality p6[ . 

The properties of (TMTSF) 2 PF6 at 10 kbar are particularly difficult to understand. For example, when the 
magnetic field is perpendicular to the current direction the magnetoresistance is smaller than when it is parallel, the 
opposite of what one observes in (TMTSF)2C104 and in conventional metals. Recently, the temperature dependence 
of the magnetoresistance when the field direction was fixed at the first magic angle was measured [Q . It was found 
to be non-monotonic: as the temperature decreased down to T m i n , the magnetoresistance decreased, it increased 
until T max , and then decreased. It has recently been proposed that these two temperatures actually represent phase 



transitions between metallic and insulating phases |26|. The magnetoresistance of the quasi- two-dimensional metal 
a-(BEDT-TTF) 2 MHg(SCN) 4 [M = K,Rb,Tl] also exhibits unusual temperature and angular dependence 

The purpose of this paper is to clarify what properties of the magic angle effects can only be explained within 
a non-Fermi liquid framework by seeing what effects can be explained within Fermi liquid theory. The interlayer 
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magnetoresistance is calculated within the framework of semi-classical transport theory. It is found that if one takes 
into account the finite bandwith along the most-conducting direction then dips in the magnetoresistance are observed 



for p/q — 1, 2, 3, ... 19 1 . Furthermore, if one assumes a simple Fermi liquid form for the temperature dependence of 
the scattering rate then at the magic angles the interlayer magnetoresistance does have a non-monotonic temperature 
dependence. Hence, one should be cautious about associating maxima and minima in the temperature dependence 
with metal-insulator transitions. However, the results obtained gives a poor description of the observed properties 
when the field is close to the b-axis and of the resistivity within the layers. 

II. CALCULATION OF THE INTERLAYER CONDUCTIVITY 
A. Semi-classical transport theory 

If the electronic dispersion relation is e(k) then the electronic group velocity perpendicular to the layers is v z = 

k dg?fc fe ^ • The interlayer conductivity can be calculated by solving the Boltzmann equation in the relaxation time 
approximation leading to Chambers' formula p{J] 

<r„ = [ v z (k)v z (k) (-^) d 3 k , (2) 
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where /(e) is the Fermi function and r is the scattering time which is assumed to be the same at all points on the 
Fermi surface. v z (k) is the electron velocity averaged over its trajectories on the Fermi surface, 

v z {k) = ~J cxp (J-^j v z (k(t))dt (3) 

where k(0) — k. The time dependence of the wave vector k(t) is found by integrating the semi-classical equation of 
motion 



dk e 



= —a V fc e x B , (4) 



If the temperature is sufficiently low that T « Ep then in Equation (Q) can be replaced by a delta function at 
the Fermi energy and Equation (||) becomes 

°z2 - f v z (k)v z (k)S(E F - e(k))d 3 k . (5) 



B. Dispersion relation along the chains 

In the tight binding approximation the dispersion relation in an orthorhombic crystal can be written as 

e(k) = —2t a cos(k x a) — 2^cos(fc„6) — 2i c cos(fc z c) , (6) 

where t a ,tb, and t a are the inter-site hopping integrals along the different crystal axes. In the Bechgaard salts, 
t a » tb,t c , the dispersion along the chains can be linearized giving e( k ) = hvp(\k x \ —kp) — 2tbCOs(bk y ) — 2t c cos(bk z ), 
where v F — 2t a asm(akp)/h is the Fermi velocity and kp is the Fermi wave vector. This linear dispersion has been 
used in a number of papers on the magic angle effect p0|- p2| , p4] | . If one solves for the interlayer conductivity within 
semi-classical transport theory one obtains 

a° zz l + ( WcO Tsin0) 2 ' 1 1 

where tv c o = evpcB/h is the frequency at which an electron traverses the Brillouin zone in the c-direction when the 
field is parallel to the b-axis. Clearly this is a smoothly varying function of and does not exhibit any magic angle 
effects. 

We now show that if the full nonlinear dispersion (|(j) is used then one does obtain magic angle effects. We will 
re-derive a result obtained earlier by Maki [Il9l with a view to elucidating the physics in the process. 
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C. Solution of the semi-classical equations of motion 



The group velocity for the dispersion relation (0) is 



v{k) 



2at a sin(afe x ) 
2bti, sm(bk y ) 
2ct c sin(cfc z ) 



(8) 



The rate of change of the wave vector k(t), in a magnetic field in the b-c plane, B = (0, B sin#, B cos 9), is given by 



^1 ^ ; — 2beBtb cos 9 s'm(bk y j 
— = — I 2aeBt a cos9s'm(ak x ) 
" -2aeBt a sin9 sm(ak x ) 



(9) 



where terms involving t c have been neglected. This is valid provided i c sin(9 <C tb cos 9. Hence, the results below will 
not be valid as 9 -> 90°. 

In order to calculate the z-component of the velocity one needs to obtain k z (t), which Equation (^c) shows is 
determined by k x (t). To zero-th order in tb, k x (t) = kp. Integrating Equation (^|b) then gives 



where 



k y (t) = ky(0) + , 



LUb = v peBb cos 9 /fi = L0bocos9 



(10) 



(11) 



is the frequency at which an the electron traverses the Brillouin zone in the direction of the b-axis. 
Substituting this into Equation (|9^i) and integrating gives, to first order in tb/t a , 

2tb 

k x (t) = kp + - — cos(6fc„(0) + ojbt) . 
nvF 

We obtain k z {t) by using Equation (||c) and substituting in (|T^), giving 
—2aeBt„ sin9 



dk z 
dt 



sin(kp) cos [ ^ cos(bk y (0) + oJbt)) + cos(akp) sin f^^l cos(bk y (0) + ujbt) 



(12) 



(13) 



\nvp j \hvp 

where we have used trigonometric identities to expand s'm(ak x (t)). If we take a linear dispersion relation, the second 

VF 

n 



term in (|13|) will equal zero and we are left with = Bev £ Sln 6 , where we have assumed that at t = 0, the wave 



vector in the x-direction (k x ) is equal to kp. 
Now to first order in tb/t a , 



dk z 
~~dt 

Integrating this we obtain 



-2aeBt a sin9 [ . . . . , f 2atb 
sm(akp) + cos{akp) sin I — — cos(ofcj,(0) + ujbt) 



k z {t)c — k z (0)c — u> c t — 70 tan6'sin(6fc J( (0) + u>bt) , 



(14) 



(15) 



where 



and 



lu c = uj c q sin ( 



2ct b a 
7o = ~z — T cot(afc F ) . 
nvp b 



(16) 



(17) 
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D. Evaluation of the interlayer conductivity 



Substitution of ( |15| ) into the z-component of the velocity gives 

v z (k z (0), (/>,</>') = sin ( ck z (0) + —<p - 70 tan#sin( 
n \ LL>b 

where <p = —LUbt, (f> = bk y (0). 

The interlayer conductivity given by (|J) can then be written in the form 



) » 



(18) 



A7T 3 b?lVF 



7r/c p27T 

dk z {0) / dcbv z (k z (0),< 
-7r/c Jo 



d<j> 




■ exp 



TUJb 



v z (k z {0), </>,<!> ) 



(19) 



We now expand Equation ( |lq ) using trigonometric identies and substitute the Bessel generating functions to obtain 

00 

Jo (70 tan0) + 2 -hk(lo tan0) cos((2/c)(0 - </>')) 
fc=i 

00 

2^ J 2 fe+i(7otan0)sin((2/c + l)(0-0')) 



sin cfc z (0) H 7 



cos I ck z (0) 



fc=0 



(20) 



and w z (fc z (0), (/>) is obtained by setting 0=0. Substituting these expressions into Equation ( fl9| ) and performing the 
integrals over (f> , (j) and dk z (0) the final expression for the conductivity becomes 



J (7 tang) 2 
1 + (w c r) 2 



^ J,, (70 tan ( 



1 + T 2 (lU c — UJbv) 2 1 + T 2 (UJ C + UJbVy 



(21) 



where a zz = is the interlayer conductivity in zero field. Note that for fixed v and z <C 1 



Ju(z) 



i> + i)' 



(22) 



Maki H obtained a similar result, although he included the corrections to « h and Ue to next order in . This 

raises the general question of to what order in tb/t a is the above expression for a zz valid. We only calculated k z (t) to 
first order in tb/t a . Strictly speaking, this means that (|l]) is valid to second order in tb/t a . However, we anticipate 
that a general solution for v z (t) will be of the form v z (t) ~ sin(w c i) J2 n a « sm ( w bi) where a„ is of order {tb/t a ) n . This 
means that the coefficients in ( f2~i| ) for v > 2 will change but be of the same order. 



III. MAGIC ANGLES 

The angular dependence of the interlayer resistivity given by Equation ( pl| ) is shown in Figure [l] for several parameter 
values. Dips occur at the "magic angles" given by io c = vuJb or 

tan6»=-^ ±i/ = l,2,3,.„. (23) 

c 

where b and c are lattice constants. The size of the dip at the v-th magic angle, compared the background magne- 
toresistance will be of order 

7o b \ 2u fuj c r\ 2 



2 /) Or) ■ < 24 > 

The size of the dips is determined by the parameter 70, defined by (p"7|), which is determined by the geometry of the 
Fermi surface. Note that if 70 — > 0, ( ^l| ) reduces to (0). This is because the limit 70 — > corresponds to taking a 
linear dispersion relation. If 70 1 then the dips will decrease in magnitude rapidly with increasing v. For example 
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if 'yg ~ 0.1 then the v = 1 feature will be five orders of magnitude smaller than the v = 3 feature. Note that when v 
becomes sufficiently large this will no longer be valid because 70 tan 9 ~ 1. 

We now consider what is a realistic value for 70 for the (TMTSF) 2 X materials. If we look at the the form of 70 in 
Equation ( p^) we note that the factor |^ equals the parameter 7 which determines the periodicity of the Danner 
oscillations^,^. For (TMTSF) 2 C10 4 if was estimated to be 0.24. The lattice constants for (TMTSF) 2 PF 6 are 
a = 7.3A, b = 1.1 A and c = 13.51, while for (TMTSF) 2 C10 4 , b is twice as large due to anion ordering [[§. The 
cot(afcir) term depends on the band filling. At three-quarter filling kp = §£, and cot(afcp) = 1< This gives a value for 
'Jo(PFq) = 0.24 and 7o(C7C>4) = 0.12. Note that for half-filling kjr = and thus 70 = and there will be no magic 
angle effects unless we solve the semi-classical equations to higher order in tf,. 

Figure [l] is qualitatively similar to experimental results for (TMTSF)2C104 at ambient pressure Q and at 6.0 
kbar |7| and (TMTSF) 2 PF 6 at 6.0 kbar (0.3 K and 4 T) Q. A small difference is that the experimental data 
shows a small dip near 90°, whereas the theoretical curve shows no such dip. It is quite possible that the small dips 
can be explained within semi-classical transport theory if one includes the effect of a finite t c in the solution of the 
semi-classical transport equations. An analogous effect occurs when the field is rotated in the a— c plane: for coherent 
interlayer transport with finite t c a peak in the angular-dependent magnetoresistance occurs when the field is parallel 
to the a axis |3l[| . 

The angular dependence of the interlayer magnetoresistance of (TMTSF) 2 PF 6 at pressures of about 10 kbar is quite 
different from that shown in Figure [j] At fields less than one tesla the angular dependence is similar to that given by 
Eqn. (0). However, above one tesla, p zz ~ (i? cos 0) 1-3 and so a large dip is observed near 90° and at 90° the in-field 
resistance is comparable to the zero-field resistance . 

A number of theoretical papers ]ll]j20| |22| have predicted effects when tan 9 = ~ where p/q is fraction. In 
contrast, the model considered here only gives effects for q = 1. A review of the experimental literature shows that 
the only reproducible fractional features seen have been in (TMTSF) 2 C104 at p/q — 3/2 and 5/2 |l2). This can be 
explained within the framework considered here. If (TMTSF) 2 C104 is slowly cooled anion ordering occurs and the 
lattice constant in the b direction doubles so in ( |23| ) b should be replaced by 2b. However, if a sample which is not 
completely anion ordered it will produce features at angles corresponding to half-integers for a fully anion-ordered 
sample. 

IV. TEMPERATURE DEPENDENCE OF THE INTERLAYER MAGNETORESISTANCE AT THE 

MAGIC ANGLES 

Suppose that the field direction is fixed at a magic angle and the temperature (and thus the scattering time r) is 
varied. 

The first magic angle {v — 1). Setting lu c = uj-^ and using the fact that 7 C 1 to take just the first term in the 
series, i.e. v = 1, the conductivity is 



a zz (9 1 ) ~ At 



2 

1 / 70 tan 9\ 



1 + (lu c q sin 9\t) 



(25) 



where 9\ represents the first magic angle and a zz — At, where A is the ratio of the zero field conductivity to r. A 
plot of interlayer resistivity verses -j^ is shown in Figure || for different values of ujo. The interlayer resistivity is a 
non-monotonic function of r. It will be seen below that this leads to non-monotonic temperature dependence. 

We now find for what values of r the maxima and minima seen in Figure ^ occur. Finding the extrema of 
Equation ( p5| ) as a function of r gives that the minimun occurs when 

^coT ^ — ^- (26) 

and the maximum occurs when 

2 

^coT- T—z-- (27) 

7o sin 9 X 

If 70 -C 1, then the maximum will only be observed at sufficiently high fields and in high purity samples. 

The second magic angle (y = 2). To obtain the interlayer conductivity for the second magic angle then a similar 
argument to that given above leads to 
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a zz (6 2 ) ^ At 



1 + (tuj c ) 2 



70 tan 9 2 



(7 tan6> 2 ) 
64 



(28) 



where we have set uj c = 2u>b and # 2 is the v — 2 magic angle. For small 70, the minima is again given by lo c t ~ 1. To 
find the maxima we expand the first two terms in ( |28| ) to fourth order in ^ r J p . The maxima occurs when 



sin 6*2(70 tan 6> 2 ) 2 ' 



(29) 



Since this is smaller than ( p7| ) by a factor of I/70, in order to see this maximum even higher fields and lower 
temperatures will be required than for the maximum associated with the first magic angle. 

Conductivity as 8 — > 90°: We can expand the term ( 1 
second order in v cos 9 to obtain 



1 + cr 



t~, — n — 2? 1 ) in the summation in (Ell) to 

(30) 



■{l + v 2 cos6 2 A + ) 



a 2 (3a 2 1) 

where a = luqt and A = ~ 3 for ljo'T" ^> 1- Substitution of this into the conductivity gives 



<r„(0->9O°) 



1 



(w r) 2 

This can be simplified using the identities 



J (70 tan 6) 2 + ^(70 tan 6») 2 (1 + 3^ 2 cos6> 2 



E ^V„(z) 2 = z 2 /2 



(31) 



(32) 



to give 



Pzz(8 = 90°) (lo c0 t) 2 



rzz 



l + 37o 2 



(33) 



The resistivity is quadratic in field as for the case of a linear dispersion, but the co-efficient is smaller. This is consistent 
with Figure [l] which shows that the resistivity near 90° does decrease with increasing 70 . This field dependence is 
quite different to what is observed in the (TMTSF) 2 X materials. For (TMTSF) 2 C10 4 at 6.0 kbar a linear field 
dependence is observed at high fields. In (TMTSF) 2 PFg at pressures from 6 to 10 kbar the resistivity saturates as the 
field increases 0,|8|. However, caution is in order because derivation of the quadratic dependence involved assuming 
that t c tan 9 <C i& and so we only expect a quadratic dependence slightly away from 90° or in the limit t c — > 0. 

Fermi liquid model for the temperature dependence: We now consider a specific model for the temperature depen- 
dence of the scattering time r. In a Fermi liquid the scattering rate, at temperatures much less than the Fermi 
temperature, has a temperature dependence of the form p6|,p2| 



to 



+ PT 2 , 



(34) 



where the first term is due to impurity scattering and the second is due to electron-electron scattering. Using this 
expression for r in Equation ( p5| ) we can now plot the temperature dependence of the resistivity. This is shown in 
Figure [| for various values of uj c qTq . The resistivity is not a monotonic function of temperature but has a minimum 
when lo c t(T) ~ 1. If the field is sufficiently high there is also a maximum. Using (34) we see that the minimum occurs 
at a temperature 



Lo c o sin 9 



1/2 



If T fnax > VO^o) then (E7|) implies 



Tn 



w c o7o sin 9 tan t 



1/2 



(35) 



(36) 
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We are unaware of any measurements of the temperature dependence of the interlayer magnetoresistance of 
(TMTSF)2C104 at the magic angles. Although the temperature dependence at the v = 1 magic angle shown in 
Figure|| is similar to that reported in Ref. jjj for the intralayer resistance of (TMTSF^PFg at 9 kbar, the observed 
temperature dependence of the interlayer magnetoresistance is different [ [33| . It depends weakly on the temperature 
from 15 K down to about 3 K and then decreases. 

The temperature dependence shown in Figure || for 9 = 90° i.e., as the magnetic field is aligned with the b axis 
the temperature dependence is qualitatively similar to that observed in (TMTSF)2C104 at ambient pressure p3. 
Again, qualitatively very different behavior was observed § in (TMTSF) 2 PF 6 at 10 kbar. There it was found that 
the in-field resistance had a temperature dependence similar to the zero-field resistance. 

Zheleznyak and Yakovenko have given a heuristic argument as to the origin of the maximum and minimum tempera- 
tures seen in the intralayer resistance in Ref. , suggesting that metal- insulator and insulator- metal phase transitions 
occur as the temperature passes through these values p6j. They argue that T min w uj\, ~ B and T rnax w t c which 
is independent of field. Above T m i n the system is a two-dimensional metal. Below T m i n , a magnetic field causes the 
electron motion in the b direction to be quantized resulting in a one-dimensional dispersion and correlations produc- 
ing insulating behavior. Below T max , the interlayer coupling becomes important and metallic behavior is recovered. 
In contrast we find that T m i n and T max are given by ( j35| ) and (|36|), respectively, and both scale with \[~B. Careful 
measurements should be able to distinguish between these two different field dependences. 



V. CONCLUSIONS 



This paper only considers the interlayer resistivity p zz , whereas magic angle effects are also seen experimentally in 
the intralayer resistivity p xx . Maki pointed out that Jl9| the semi-classical theory will only give resonances in p xx of 
order (t c /t a ) 2 whereas they are observed to be much larger. A possible way around this problem is that experiments 
that are meant to measure p xx may actually be measuring some of p zz . This is because in highly anisotropic metals 
it is difficult to arrange the contacts and current path so it lies completely within the layers. This potential problem 
increases the motivation to make thin film samples of these metals. 

It has been shown that within semi-classical transport theory a nonlinear dispersion parallel to the chains is 
necessary to produce dips in the interlayer magnetoresistance at integer magic angles. If the field direction is fixed at 
one of the magic angles then one observes both minima and maxima in the temperature dependence of the interlayer 
magnetoresistance. This arises from the temperature dependence of the scattering rate. Since these maxima and 
minima can exist within a Fermi liquid model one should be cautious about associating them with non- Fermi liquid 
behaviour or metal-insulator transitions. On the other hand, although the Fermi liquid model considered here gives 
a good description of many of the properties of (TMTSF^PFg at pressures from 6 to 8 kbar and of (TMTSF^CICU 
it gives a poor description of their properties when the field is parallel to the layers and of the intralayer transport. 
Qualitatively very different behavior is observed in (TMTSF^PFg at pressures of about 10 kbar; explaining it remains 
a considerable theoretical challenge. 
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FIG. 1. Angular dependence of the interlayer magnetoresistance of a quasi-one-dimensional Fermi liquid. 9 is the angle 
between the magnetic field and c-axis and the magnetic field is rotated in the b — c plane. The dips in the resistivity occur at 
the magic angles defined by tan 6 = v where, v = 1 and 2 (see Equation (p3|)). The dips only occur when one takes into account 
the non- linear dispersion parallel to the chains and their intensity depends on 70, which is determined by the geometry of the 
Fermi surface (Eq ( |l7| ) ) . r is the scattering time and ujq is the frequency at which the electron traverses the Fermi surface when 
the field is perpendicular to the layers. The lattice constants 6 and c are set equal and so u>o = cuto — uj c q. The resisitivity is 
normalized to p° zz , the interlayer resistivity at zero magnetic field. 




FIG. 2. Non-monotonic dependence of the interlayer resistivity, at the v = 1 magic angle, on the scattering time r. The 
curves shown are for 70 = 0.25 and for various values of ujq, which is proportional to the magnetic field. 1/s/t is used for the 
horizontal axis because it will an increasing function of temperature in a Fermi liquid. A local minima occurs at ljqt ~ Bi ^ g 
and there is a local maxima when loot ~ 70 3 ? n e ■ The interlayer resistivity p zz is normalised to -jj, where A is a constant equal 
to the ratio of the zero field conductivity to the scattering time r. 
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FIG. 3. Non-monotonic dependence of the interlayer resistivity, at the v = 1 magic angle, on temperature. A Fermi liquid 
form for the temperature dependence of the scattering rate is assumed. A value of to/3 = 0.025 K~ 2 is used in Equation (134) 
so that the temperature dependence of the zero-field resistivity roughly corresponds to that of typical samples of Bechgaard 
salts. At a temperature T min defined by wo T (T m in) — 1 there is a minimum in the resistivity. For sufficiently high fields there 
is a temperatute T ma x at which the resistivity is a maximum. 
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